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Abstract 

A mean field model is used to investigate if a plateau in caloric curve can 

be reached in mononuclear configuration. In ttie model the configuration will 

break up into many pieces as the plateau is approached. 
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The appearence of a plateau in the experimental nuclear caloric curve, since it was 
pointed out first in [1] , has continued to be an important issue in intermediate energy heavy 
ion collisions [2]. The plateau would signify a maximum in specific heat and could be a 
signature of a phase transition. 

The appearence of a maximum in the specific heat at temperature T ^ 5 MeV was seen in 
theoretical calculations much earlier [3]. The model is the SMM (Statistical Multifragmen- 
tation Model) where it is assumed that the nucleus breaks up into many pieces at a volume 
significantly larger than the normal nuclear volume. Subsequent modelling has reinforced 
this picture of a heated nucleus breaking up into pieces as the phenomenon which gives a 
peak in the specific heat. Take, for example, the LGM (Lattice Gas Model) [4]. Here at low 
temperature there is a large percolating cluster. As the temperature is raised one reaches 
a point when the percolating cluster breaks up into many smaller clusters. It is here that 
the maximum in specific heat is seen (Fig. (19. 2) in [2]). This is more dramatically demon- 
strated in the canoncal thermodynamic model (which is in the same spirit as the SMM but 
is much easier to implement). Here also at low temperature there is a large blob of matter 
which breaks up into many pieces and again the maximum in specific heat is obtained in 
this region. The thermodynamic model can be extrapolated to large numbers of particles. 
In this limit it is shown that the break up is very sudden as a function of temperature. This 
is a case of first order phase transition and the maximum of specific heat is obtained at the 
phase transition. 

The present work is inspired by a recent calculation which showed that, in a different 
model, a plateau in the caloric curve (hence a maximum in specific heat) can be reached 
in mononuclear configurations as well [5]. One assumes that e*(=the excitation energy per 
particle) as a function of density increases quadratically about the ground state density. 
At a given excitation energy, the nucleus expands in a self-similar fashion till it reaches its 
maximum entropy. Effects of interaction on entropy is taken through a parametrisation of 
m * /m where m* is the effective mass. A temperature is defined microcanonically. The 
authors then find that when they plot temperature against excitation energy a plateau is 
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found around 5 to 6 MeV. 

A more familiar model in nuclear physics which allows study of the caloric curve in 
mononuclear configurations is the temparature dependent mean field model (Hartee-Fock 
and/or Thomas- Fermi model). This has certain advantages. When one does a standard 
mean field calculation at a fixed temparture, one minimises the free energy F — E — TS 
[6]. This means that when we get the self- consistent solution at a given temperature, we 
have obtained a solution which has zero pressure. If we draw a caloric curve with energies 
of these solutions this caloric curve pertains to zero pressure. The specific heat that we will 
get will be Cp with p = 0. 

Investigation of the caloric curve with temperature-dependent Thomas- Fermi theory was 
done in the past [7]. For nuclei ^^'^Sm and ^^Kr caloric curves were drawn and a maximum in 
specific heat at temperature 10 MeV was found. But there is ambiguity whether the systems 
stay mononuclear or not. Even though we are interested in finite nuclei, when we do finite 
temperature mean field theory beyond a certain temperature, fermi occupation factors for 
orbitals in the continuum will cause the finite system to spread out. The calculations have 
to be done in a big box and the size of the box will affect the answers. Although given 
the box, definite numerical results can be obtained it will be hard to decide whether it is a 
mononucleus or a system of gas. This is best explained using figs. 1 and 2 of [7]. Looking 
at figures one would venture to say that at temperature 5 MeV we have a mononucleus but 
at temperature 10 MeV the nucleus spreads to 12.5 fm with density «i 0.008/m~^. Is this a 
gas or a nuclear liquid? 

A much more definite answer can be obtained in the nuclear matter limit. This is 
pursued in this work. For a fixed temperature we do calculations for different densities. The 
density where the free energy per particle is minimised is the solution for this temperature; 
e* of this solution is the appropriate e* for this temperature. As expected, starting from 
zero temperature, the system expands. The minima of free energy drop to lower and lower 
density as the temperature increases. But beyond a certain temperature, minimum in free 
energy disappears. The nucleus will now break apart. This happens before T flattens out 

3 



as a function of e*. This is shown in figs. 1 and 2. 

It remains to give some details of the calculation. We use the mometum dependent mean 
field of [8,9]. The potential energy density is given by 

= ±t + + £ / jaw'^jmm (1) 

2a, <t+1 p5 p„J J ^ ^ i + l'^P 

Here f{f,p) is phase-space density. In nuclear matter, at zero temperature f{f,p) = 
■^Q{pf — p) where 4 takes care of spin-isospin degeneracy. At finite temperature the theta 
function Q(pF—p) is replaced by Fermi occupation factor (see details below). The potential 
felt by a particle is 

Po Po PqJ 1 + [^]^ 
Here A =-110.44 Mev, 5=140.9 MeV, C=-64.95 MeV, po = O.lG/m-^, a=l.2A and A = 
1.58p^. This gives in nuclear matter binding energy per particle=16 MeV, saturation density 
Po — .16 fm"^, compressibility K—215 MeV and m * /m=.67 at the fermi energy; u{p,p) 
gives the correct general behaviour of the real part of the optical potential as a function 
of incident energy. A comparison of u{p,p) with that derived from UV14-I-UVII potential 
in cold nuclear matter can be found in [9]. The specific functional form of the momentum 
dependent part arises from the Fock term of an Yukawa potential. Mean fields given by eqs. 
(1) and (2) have been widely tested for fiow data [10] and give very good agreement. 

To do a finite temperature calculation the following steps have to be executed. We need 
to find the occupation probability 

= e^Krt-Ml + 1 

for a given temperature 1/(3 and density p. If e(p) were known a priori, this would merely 
entail finding the chemical potential from 

But the expression for e(p) is 



<P) = ^ + A[^] + B[^r + RM (5) 
2m po Po 

where at finite temperature 

Thus knowing R{p,p) requires knowing n[e(p')] aheady for all values of p'. This self- 
consistency condition can be fulfilled by an iterative procedure (details can be found in 
[9]). 

To calculate pressure, we use the thermodynamic identity pV — —E + TS + jiN which 
then gives 

p = a + h + c (7) 

Here 



167r i-cc, 

a = —V — 



where v is given by eq. (1) and the second term is the contribution from the kinetic energy. 

IGtt r'^ 

b=—T—— / p^[nlnn + (1 — n) ln(l — (9) 

The term c is up. The free energy per particle is — (a + b)/p. The test p = when the 
minimum of free energy is reached provides a sensible test of numerical accuracy. 

A second set of calculation was done with a potential energy density without any momen- 
tum dependence. This means the C terms in eqs. (1) and (2) are put to zero and A, B and 
(7 are readjusted to give desired values of binding energy per nucleon, equilibrium density 
and compressibility. The constants now are A=-356.8 MeV, 5=303.9 MeV and a ==7/6. We 
study equilibrium values of density with temperature with both the interactions but so, as 
not to clutter up figures, only the curves with momentum dependent interaction is shown 
in fig.l. In fig. 2 we draw the caloric curves with both sets of interaction and for reference, 
Fermi gas results e* = aT^ are also shown. 



In [5] , the effects of interactions on the caloric curve were taken through an effective mass 
m* replacing the real nucleon mass m. Normally one writes m* /m = rrikm^^. In momentum 
independent mean field both mu and are ignored and m * /m=l. In the momentum 
dependent mean field calculation that we have done here is in but is not. But in 
both the models, the nucleus will break up before T flattens out against e*. As in [5], with 
momentum dependence the curve is initially above the Fermi gas and meets the Fermi gas 
curve at some later point. 

The mean field model thus reinforces what has long been perceived in other models. At 
low temperature, the system has a large mass. There are probably some lighter masses too. 
As the temperature increases the larger mass which has the bulk of the system will break 
up. This is the point where the caloric curve flattens out. We also note that experimentally 
the interesting features of the caloric curve occur around 6 MeV temperature whereas in figs 
1 and 2 they are occurring around 12 MeV temperature. The finiteness of the system will 
bring this down as will the Coulomb interaction. In any case some overestimation in mean 
field theory is expected. In the Ising Model mean field calculations overestimate the critical 
temperature by fifty per cent [12]. 

Lastly the lessons from the simple mean field model done here can be used to bolster 
the expanding emitting source (EES) model of W. Friedman [13]. The model assumes that 
initially the hot system evaporates as well as expands. For low initial temperature the 
system will cease expanding and will revert towards normal density. But beyond a certain 
temperature at the end of this slow expansion the system will explode. 

I thank Lee Sobotka for many communications. This work is supported in part by the 
Natural Sciences and Engineering Research Council of Canada and in part by the Quebec 
Department of Education. 
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FIG. 1. The free energy per particle as function of density at different temperatures. This is 
plotted as a difference from the minimum value of F in the frame. At temperature 12 MeV there 
is no minimum and the system will roll down to lower density. For this temperature already at 
lower values of p in the figure, the derivative ^ is negative, i.e., the system has entered a region 
of mechanical instability. 
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FIG. 2. Temperature against e* in mean field model with momentum dependence (MFl), in 
Fermi gas model (FG) and in mean field model without momentum dependence (MF2). As in ref. 
5, the curve with momentum dependence lies above the Fermi gas curve till they meet (around 
temperature 11.5 MeV). The MFl curve has not yet flattened out. Just below temperature 12 
MeV the curve stops. The nucleus will break up into many pieces at higher temperature. Similar 
situation happens with MF2 except that here it will break up many pieces at ~ 12.5 temperature 
again before the caloric curve flattens out. 
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